Mixings of 4-quark components in light non-singlet scalar mesons in QCD sum rules 
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Mixings of 4-quark components in the non-singlet scalar mesons are studied in the QCD sum rules. 
We propose a formulation to evaluate the cross correlators of qq and qqqq operators and to define 
the mixings of different Fock states in the sum rule. It is applied to the non-singlet scalar mesons, 
ao and Kg. It is found that the 4-quark operators predict lower masses than the qq operators and 
that the 4-quark states occupy about 70-90% of the lowest mass states. 

PACS numbers: 



INTRODUCTION 



Hadrons which are not represented by minimal qq (me- 
son) or qqq (baryon) configurations are called exotics, and 
+ with baryon number B = 1 and strangeness S = +1 
is the most prominent example. Recent developments in 
hadron spectroscopy are largely enhanced by the discov- 
ery of the pentaquark + [lj. The LEPS group, who first 
observed a peak in the NK spectrum produced by the 
SPring-8 photon beam, has confirmed the resonance in 
their succeeding experiment 0, while many others have 
been unsuccessful in producing the resonance in various 
types of experiments [3j . While experimental efforts are 
obviously needed now, a large number of theoretical stud- 
ies have been carried out to explain the mass and the 
width of + , as well as its structure and reactions. At 
the same time, various new types of exotic hadrons have 
been proposed and examined. In fact, the spectroscopy of 
old (ordinary) hadrons has been reexamined, and various 
possibilities of multiquark components of hadrons have 
been pointed out. Among them, the possible 4-quark 
structures of the scalar mesons are not new, but rather an 
old idea 0, [5[ . The (non-relativistic) quark model based 
on SU(6) x 0(3) symmetry does not easily account for 
the mass pattern of the lowest-lying scalar-meson nonet, 
(cr(600),a (980),/ (980), J ft: o *(800)). (if * has been indi- 
cated in Ktt final states in J/tp and D-meson decays, 
but has not yet been established Q.) The expected mass- 
ordering as qq states would be m{a) ~ m(ao) < m(fo), if 



we assume ideal mixing, i.e., a 



uu+dd 



~V^\ a ° 7T 

fo ~ ss. This pattern does not agree with the one from 

experiment. Furthermore, while the scalar mesons are 

classified as 3 Pq states, their spin-orbit partners J = 1 

and 2 states are not observed in their vicinity. 



i — dd 



and 
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A possible solution to this puzzle is to consider 4-quark 
exotic states for the scalar mesons Q. Suppose that di- 
quarks with flavor 3, color 3 and spin 0, i.e., 



U = (ds) S= 0,C=3J=3, 
S = (ud) S =0,C=3,f=3, 



D = {su) S =0,C=3,f=3, 



(1) 



are the building blocks of the scalar mesons. Then the 
scalar nonet in the ideal mixing may appear as 

er ~ SS ~ (ud)(ud), 

a ~ -t=(UU ~ DD) - -j= [(ds)(ds) - (su)(su)] , 

V ^ 



V2 



fo ~ -L(UU + DD) ~ -L [(ds)(ds) + (su)(su)] . 

Then one sees that the strange-quark counting suc- 
cessfully predicts the observed mass pattern: m(a) < 
ro(ao) ~ Tn(fo). It also explains why the J = state 
is isolated from J = 1 and 2 partners and thus helps to 
explain anomalies in the scalar-meson nonets. There are 
other hadrons which are pointed out to have possible ex- 
otic multi-quark components, such as D*, A(3872) and 
A(1405), where the last one may be a baryon resonance 
composed of 5 quarks @, H, H, EH H3 • 

Does the QCD dynamics allow such states? We ac- 
tually have a simple reasoning why the scalar mesons 
have significant 4-quark components. As qq states, the 
scalar mesons, which have positive parity, must contain 
orbital excitation L — 1 (P-wave mesons). It is gen- 
erally accepted that the orbital excitation in the quark 
model requires additional 500 MeV, for instance, as is 
seen from the mass differences between the positive- and 
the negative-parity baryons; iV(1535) — iV(94Q). On the 
other hand, the 4-quark components, qqqq, may be real- 
ized in L = without orbital excitation and thus have 
an advantage over the P-wave excited states. In fact, the 
cost of having an extra qq pair may be about 500 MeV; 
the average mass of 7r and p mesons. In terms of the di- 
quarks, the combination of two scalar diquarks in L = 
is quite preferable in the standard quark-model and may 
easily reproduce the observed spectrum of the scalar 
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mesons. Therefore it is quite interesting and important 
to answer the question whether QCD indeed gives multi- 
quark states a lower energy and induces mixings of the 
different Fock components. 

Our purpose in this paper is to study the possible mix- 
ing of different Fock components in the QCD-sum-rule 
approach. It is important to establish a well-controlled 
formulation of the mixings of two (or more) types of op- 
erators in producing hadrons from the vacuum. We also 
examine the definition of the Fock-state-mixing param- 
eters so that we can compare the results with the pre- 
dictions of the quark models. In this paper, we focus on 
the 4-quark components of the flavor non-singlet scalar 
mesons, ao and Kq, and demonstrate how our formula- 
tion works in this case. 

In Sect.ini we present the basic ideas of the QCD-sum- 
rule approach [lj, EH and give the interpolating field op- 
erators corresponding to the qq and qqqq components of 
the scalar mesons. In Sect. IIIH we consider the definition 
of the mixing parameters of the two Fock components. 
Two distinct methods of defining the mixing parameters 
are presented. In Sect. IIV1 we give the results of the 
sum rules and study their significances. The conclusion 
is given in Sect. El 



II. THE QCD SUM RULES 

The sum rule is obtained by expressing the two-point 
function 

n(p 2 ) = ijd 4 x e lpx (0\T[J(x) J f (0)]|0), (2) 

in two ways. One of them is based on the operator prod- 
uct expansion (OPE), where Eq. © is calculated in deep 
Euclidean region (— p 2 — > oo) and is described in terms 
of the QCD parameters, such as the quark condensate 
(qq), the gluon condensate (a s ir~ 1 G 2 ) , the current-quark 
masses m q and so on. The other one is based on a phe- 
nomenological parametrization of the spectral function. 
The spectral function at the physical region (p 2 > 0) is 
assumed to have a sharp-peak resonance at p 2 = m 2 and 
continuum at p 2 > Sth- 

Pphcn.b 2 ) = \X\ 2 S( P 2 - m 2 ) + 0(p 2 - s th )poMp 2 )- (3) 

The sum rule is obtained by matching the two expressions 
so that the mass of the resonance, m, and the other phe- 
nomenological parameters can be determined from the 
QCD parameters. The two expressions of the correlators 
are connected by a dispersion relation. To simplify the 
sum rule, we approximate the continuum spectrum as the 
same form of the spectral function on the OPE side. 
We further apply the Borel transformation, which is, 

-p 2 ,n 2 ^oo TV. \dp z J 



where the limit is taken with M 2 = fixed. It sup- 
presses large-p 2 region by the factor of e~ p2 / M2 , and 
thus the pole contribution is enhanced. Due to the p 2 
derivatives, the subtraction terms vanish. The Borcl 
transformation also suppresses the effects of the higher- 
dimensional terms in the OPE. 

We employ the following local operators for Oq : 

J 2 (x) = (u a (x)d a (x)) , 

Ji(z) = tabctdec (dl(x)C-f 5 S b (x)) (s d (x)j 5 Cu^ (x)) , 

(5) 

where a,6, . . ., represent the color and C = ij 2 j°. The 
2-quark operator, J2, is uniquely determined, while a 
specific 4-quark operator, J4, is taken so as to consist 
of two + -diquark operators. The flavor structure of J 4 
is arranged to be UD, and thus identical to J 2- Simi- 
larly, the 2-quark and the 4-quark operators for Kq are 
given by the SUf(3) rotation from the fields for ao, i.e., 
u — ► d, d — ► s, s — ► u in Eq. (|5|). 

The number of independent 4-quark operators is dis- 
cussed in a previous paper [l4l |. There are two kinds of 
flavor-octet 4-quark operators: i.e., 3 C x 3 C and 6 C x 6 C . 
Each kind allows five independent operators; thus alto- 
gether there are ten possible operators. 

It is, however, difficult to find the most suitable op- 
erator for the QCD sum rule of the scalar mesons. In 
Ref . [lij . the authors attempt to find the best combina- 
tion. It is interesting to employ such an operator for the 
present method, while the computation may be much 
more complicated. 

We here employ, instead, the simplest one composed 
of the 3/ scalar diquarks, J p = + and C = 3, which are 
supposed to be bound most strongly. 

Our results of the OPE side are listed in Appendix [A"l 
The OPE is truncated after the dimension-6 terms for 
Il22(p 2 ) (see Eq. (|AI[) '). In order to deal with the same 
order of power corrections, we expand Ii-2i(p 2 ) (Eq. (|A2|) ) 
up to dim-9 and n 44 (p 2 ) (Eq. J55J)) up to dim-12. High 
dimensional operators, such as (qqqq) and (qqqqG) are 
evaluated by the vacuum-saturation approximation. 

We consider the diagrams containing the qq annihila- 
tion in order to take into account the mixing of different 
Fock states. The quark-pair annihilations are substituted 
by (qq) or (qg s a-Gq) (see FIGH]). Because the interpolat- 
ing fields (O have the normal ordering, the perturbative 
part of the qq annihilation must disappear. Because the 
OPE is represented as a polynomial in x, only the zero-th 
order term survives in the x — > limit. 



III. MIXING OF FOCK STATES 

The main goals of this study are to determine whether 
the scalar mesons are dominated by qq or qqqq, and if 
they are mixed, to calculate their probabilities. It turns 
out that such a mixing is not easily quantified. It would 
be natural to consider the strengths of the couplings of 
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FIG. 1: The annihilation diagrams. 



the 2-quark and 4-quark operators to the physical state 
and then evaluate the mixing angle. However, such a 
procedure is strongly dependent on the definition and 
normalization of the local operators. Indeed, a numerical 
factor can be easily hidden in the local operators and thus 
the magnitudes of the coupling strengths are ambiguous. 

This problem happens to be more fundamental than 
just the definition of the operators. In quantum mechan- 
ics, mixings of Fock states may occur if the Hamiltonian 
contains terms which change the number of particles. 
Then the wave function of an eigen-state can be written 
uniquely as a linear combination of the Fock components, 
each of which is normalized properly. In the field theory, 
however, the Fock-space separation may not be unique as 
the number of quarks, which is defined by the number of 
quarks + number of antiquarks, is not a conserved quan- 
tum number. One may not be able to "measure" the 
number of quarks without ambiguity because no con- 
served charge corresponding to the number of quarks 
is available. Thus we have to consider the concept of 
the "number of quarks" in the context of the quantum- 
mechanical interpretation of the field-theoretical state. 

Here, we propose two ways to define the ratio of the 
Fock-space probability by taking the mixing of qq and 
qqqq in the scalar (arj) meson as an example. 

In the first approach, we define local operators "nor- 
malized" in the context of a full 4-quark operator J 4 in 
Eq. ([5]). In fact, the 4-quark operator J4 contains qq 
component, which is obtained by contracting the qq pair 
by the quark condensate: 



Ji{x) = J' 4 (x) + -(ss)J 2 (a;) 

J!, 



(G) 



We regard J 2 and J 4 as "normalized" 2-quark and 4- 
quark fields, respectively. The quark condensate gives 
the dimension of normalization. Using J 2 (x) and J^(x), 
we define the mixing parameter, 9. by 



<o|4(z)K) 

<0|J£(aO|oo) 



= XcosO (j){x), 
= \sin0 d>(x). 



(7) 



This definition happens to be equal to defining so that 
J a (x) = cos0J' 2 {x) + sin0J' i {x) couples to the physical 
state most strongly: 



\J a (x)\a } = A0(x), 



(8) 



where <j>{x) denotes the wave function of the center-of- 
mass motion of a (i.e., a plane wave for a momentum 
eigenstate). 

Then the mixing parameter can be evaluated from the 
following three correlation functions with an assumption 
that the poles are at the same position: 

-Imi f d 4 xe lf>x (0|T[4(x)(4) t (0)]|0) 



= S(p 2 - m 2 )|A| 2 cos 2 + cont., 
-Imi / rf 4 xe ipx (0|T[4(x) (4)^0)110) 

= 5(p 2 - m 2 )|A| 2 sin 2 + cont., 



-Imi / d 4 x e lpx (Q\T[J' 2 (x){J'^ '(0)]|0) 

7T 



5(p 2 



2 )|Af 



, sin 20 



cont.. (9) 



This definition of the mixings is model independent, but 
it depends on the choice of the local operators. Therefore, 
it does not necessarily have a direct relation to the mixing 
parameters employed in the quark models. 

In order to define a mixing angle more appropriately to 
the quark models, one may determine the normalization 
of the operators using quark-model wave functions. For 
instance, the local operators in Eq. ([5]) , can be normalized 
to the wave functions of qq or qqqq states in the MIT bag 
model: 



\(m)p ) 



V2 



(\PS) + |SP» 



^(\n)-\lV)®\ud)®-^=6 ab \ab), 



\{qqqq)s ) = \ssss) 
1 



2 (I uu) - 1 tut) -um) + 1 mt)) 

■ze a/ 3 u cs p d\ai(35p) ® — =e abc e dec \abde) , 



2V3 



(10) 



where a, (3, . . . represent the flavor and a,b, . . . represent 
the color. Then we compute the matrix elements 



(0|J 2 (0)|(g<z)p ) = A 2 0(O), 
{0\,h(0)\(qqqq) So ) = A#(0), 



(11) 



where (/>(0) denotes the center-of-mass wave function of 
the bag-model state. Now, assuming the bag-model 
states (with definite number of quarks) are normalized 
properly, one can use A2 and A4 for the normalizations of 
the operators J 2 and J4, respectively. 

In calculating the mixing parameter, one needs only 
the ratio of A2 and A4, which is given by 

~._ . 24 4, 2) 

it A 



A4 i 

A 2 47r ^2(51/2)^2 (P1/2) 

ER„ 1 



Rl /2 \j (ER n )\ ^2ER n {-l + ER n )' 



4 



AA„(P 1/2 ) = 



ER n 



1 



and 



ER n = 



2.04 
3.81 



for 
for 



R n ,2 MER n )\ ^/2ER n (l + ER n y 

Sl/2 
Pl/2 ' 

where R is the bag radius and gives the dimensional scale 
of the normalizations of the two operators that have dif- 
ferent dimensions. We here assume that the bag radius 
of the qq state is same as that of qqqq state. This as- 
sumption may be necessary to consider their mixing in 
the physical state. It is, in fact, not a bad assumption 
because the bag radius of the qq state must be larger, 
because of the relative P-wave motion, while that of the 
4-quark state is large, because it has more quarks. 

Thus, the physical state for a is given by the mixings 
of the two states: 



l a o) 



icos6\{qq) Pa ) 
cos8\ao(2q)} 



f Bin 6\(qqqq) s ) 
sin0|a o (4g)). 



(13) 



Note that here the factor i is necessary to keep the phases 
of the 2- and 4-quark states in accord so that the mixing 
parameter can be defined as a real parameter. 

In the QCD sum rules, the mixing parameter with the 
bag-model normalization can be calculated as 

2 



lmfd 4 x e^(0|4(a;)(4)t(0)|0) 



tan 2 9. (14) 



Im fd^ e v x (0\J 2 (x)4(0)\0) 

This definition is model dependent, but it gives a direct 
interpretation associated with the quark model. 
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IV. RESULTS 
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FIG. 2: The masses of oo for the pure 2-quark and pure 
4-quark operators as functions of the Borel mass, M. The 
threshold parameter, y/sth, is fixed to 1.4 GeV. 



The mass formula is given by 



m 



d(-l/M») 



Bin 



OPE 



(p 2 )9(s th -p 2 )} (AP 



Bin 



OPE 



(p 2 )e(s th - P 2 )} (m 2 



(15) 



We plot the mass of clq in the case of pure 2-quark and 
pure 4-quark in FIG. [2J We use the values of the QCD 
parameters given in Table [H where (uu) = (dd) = (qq) . 
We neglect the masses of the u- and the d-quarks. We 
see that the Borel stability is fairly good. The positions 
of the poles of the 2-quark and the 4-quark are close to 
each other. The 4-quark state is slightly lighter. 
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3: Masses of ao for three choices of the threshold param- 
y/sth = 1.3, 1.4 and 1.5 GeV. The 2-quark and 4-quark 
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ctcr. 

operators are mixed 



We plot the mass of the mixed ao for various threshold 
parameter values, Sth, in FIG. [3] The predicted mass 
extracted from the mixed operator, J a (x), is about 0.9 ~ 
1.1 GeV, which is similar to the result from the pure-4- 
quark operator given in FIG. [U 
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FIG. 4: The mixing parameters of ao plotted as functions of 
the Borel mass, M. The threshold parameter, y^Sth, is fixed 
to 1.4 GeV. 

The mixing parameters from the first-normalization 
method and the normalized operators according to the 
bag-model wave functions are plotted for various bag 
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TABLE I: Standard values of the QCD parameters. 





{qq) 


(ss) 


ml = (qg s o ■ Gq)/(qq) 




Us 


0.12 GeV 


(-0.23 GeV) 3 


0.8 x (qq) 


0.8 GeV 2 


(0.33 GeV) 4 


0.4 



radii in FIG. [4] We observe that the mixing parame- 
ter is almost independent of the Borel mass, M. We 
find that it is independent of threshold parameter, sth, 
too. The weak M-dependence of the mixing parame- 
ter, which is dimensionless, is attributed to the cancel- 
lation of the p-dependences (M-dependences) of the 2- 
quark and 4-quark correlators. The mixing parameter 
from the first-normalization method is about 70 degrees; 
in other words, the 4-quark component occupies 90% 
of ao- In the bag model, we take the central value as 
R = 4.8 GeV -1 , which is the bag radius determined for 
the qq component of ao [l5| . The mixing parameter for 
R = 5.2 GeV -1 is almost the same as the one in the 
first-normalization method. We see that the dominant 
component of ao (about more than 70%) is the 4-quark 



state for R > 4.4GeV" 
method. 

100 



and for the first-normalization 



80 



2-quark(V^=1.4GeV) 
4-quark(V.s, h =1.4GeV) 
mixed(V^=1.4GeV) 




FIG. 5: The pole contribution defined by Eq. (|16[l for ao is 
plotted as a function of the Borel mass, M, in the cases of 
pure 2-quark components, pure 4-quark components and a 
mixing between them. The threshold parameter, ^/s t h, is 
fixed to 1.4 GeV. 



We examine how well these sum rules work, 
the pole contribution defined by 

B[u OPE (p 2 )e(s th - P 2 )] 



Bin 



OPE 



(P 2 )} 



We plot 



(16) 



in FIG. [5l The valid Borel window is taken as the re- 
gion where the pole contribution is more than 30%. This 
constraint is weaker than the one usually adopted [HI], 
but we observe that the results are not sensitive to the 
choice of this value. In FIG. it is seen that the pole 
contribution for the pure 4-quark correlator is less than 
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FIG. 6: The masses of Kg for the pure 2-quark and pure 
4-quark operators as functions of the Borel mass, M. The 
threshold parameter, ^/s t h, is fixed to 1.4 GeV. 



the one for the pure 2-quark correlator. The reason is 
that the OPE of the 4-quark correlator contains higher 
powers of p 2 and grows rapidly for large p 2 . The sum 
rule for the mixed operator shows similar results. We 
set the Borel window M < 1.2GeV. We observe form 
FIGOthat the positions of the poles of the 2-quark and 
4-quark correlators begin to separate in this region. If 
the positions of the poles are too different, then the mix- 
ing will disappear. But in this case the results for the 
mixed operator and the 4-quark operator are smoothly 
connected. We conclude that the 4-quark component is 
dominant in ao- The results of Kq for pure 2-quark and 
4-quark components arc presented in FIG. [6] As com- 
pared to the Borel behavior of ao in FIG. [H the mass 
extracted from the 2-quark correlator for Kq is heav- 
ier, while that of the 4-quark components lies at about 
the same mass as ao- The experimental value of the Kq 
mass has been reported as about 0.84 GeV, and it seems 
to have a large width: about 600 MeV |6|. The scalar- 
nonet spectrum from our sum rule with the 4-quark com- 
ponent shows better agreement with the observed values 
than that with the 2-quark operator. The position of 
the pole of Kq for the 2-quark operators starts splitting 
at low M values from that of the 4-quark operators. We, 
however, attempt to evaluate the mixing parameter by 
assuming that the positions of the poles are at the same 
position. The results are shown in FIG. [7] and FIG. 
which correspond to FIG. [3] and FIG. 0] in the case of 
ao, respectively. One sees that the 4-quark component 
is dominant, occupying about 70% ~ 90% of Kq. These 
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FIG. 7: Masses of Kq for three choices of the threshold param- 
eter, ^/sth = 1.3, 1.4 and 1.5 GeV. The 2-quark and 4-quark 
operators are mixed. 
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FIG. 8: The mixing parameters of Kq plotted as functions of 
the Borel mass, M. The threshold parameter, %/sth, is fixed 
to 1.4 GeV. 



results for Kq are similar to those for a®. 

The pole contribution for Kq is shown in FIG. [9l We 
observe from FIG. [5] and FIG. [9] that the behavior is 
similar to that for clq. 



V. CONCLUSION 

A formulation is proposed to take into account the 
mixing between different Fock states in the QCD sum 
rules. In order to quantify the mixing probability, one 
needs the normalization of the operators. We suggest 
two ways: One is to define the "normalization" using a 
multiquark operator which couples to the ground state 
most strongly. The other way is to adjust the opera- 
tors to the normalized wave functions from the MIT bag 



FIG. 9: The pole contribution defined by Eq. lfTo) for Kq is 
plotted as a function of the Borel mass, M, in the cases of pure 
2-quark components, pure 4-quark components and a mixing 
between them. The threshold parameter, y / 3th, is fixed to 
1.4 GeV. 



model. 

We apply the formulation to ao and Kq which are 
members of the scalar nonet. Our sum rules indicate 
that 70 ~ 90% of these scalar mesons is composed of the 
4-quark components. We find that Kq is almost degen- 
erate with an.. 

There exist several other studies of the 4- quar k scalar 
nonet in the QCD sum rule. Two of them [l4l are 
consistent with our work in reproducing the lighter 4- 
quark mass than 2-quark mass. In these analyses, the Kq 
is predicted to have a smaller mass than ao- It should, 
however, be noticed that the predicted masses depend on 
the threshold parameter rather strongly. Therefore, one 
may choose the threshold so that Kq is lighter than ao 
reproducing experimental data. A'q is a broad resonance 
which is not completely established. Therefore we do 
not further examine the mass difference of ao and Kq. 
Another work [l7| claims that no signal is obtained for 
the 4-quark scalar nonet. We note that their analysis 
differs from ours in not considering the mixing diagrams, 
truncating dimension in the OPE, interpolating field and 
the definition of the coupling strength of ground state. 

It is extremely interesting that some of the excited 
hadrons are accounted for by considering exotic multi- 
quark components. There are several other "anomalous" 
hadrons which may contain exotic multi-quark compo- 
nents. Among them, the P-wave baryons, in particular 
A(1405), are strong candidates to have 5-quark compo- 
nents. In another paper p~8|, we report the results of ap- 
plying the current method of the Fock-space mixing to 
the flavor-singlet A state. Multi-quark mesons may also 
be found in heavy quark systems, where newly found 
states do not fit well in the qq spectrum and thus are 
suspected to have 4-quark components. 

It is also important to examine the decay widths and 
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the branching ratios of those multi-quark states. The 
mechanism for the fall-apart decay in which the multi- 
quark hadrons dissociate into two color-singlet hadrons 
without creating qq pairs. The widths associated with 
the fall-apart processes depend strongly on the quan- 
tum numbers as well as on the configurations of the 
multi-quarks. As QCD does not forbid such multi-quark 
states, the width is the key to understand why we do 
not see many "exotic" hadrons in nature. Their possibil- 
ity should be further pursued both experimentally and 
theoretically. 
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APPENDIX A: THE RESULT OF OPE 



The results of the OPE are summarized as: 

n 22 (p 2 ) = i J d 4 xe ipx (0\T[J 2 (x)4(0)]\0) 

3 1 - 1 

= -— -^p 2 \n(-p 2 ) - — (m u (uu) +m d (dd)) ^(m d (uu) + m u {dd)) 
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(p 2 fM-p 2 



4! 5!2 8 tt 61 
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